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Abstract We model the low-temperature specific heat of solid 4 He in the 
hexagonal closed packed structure by invoking two-level tunneling states in 
addition to the usual phonon contribution of a Debye crystal for tempera- 
tures far below the Debye temperature, T < <9_d/50. By introducing a cutoff 
energy in the two-level tunneling density of states, we can describe the excess 
specific heat observed in solid hep 4 He, as well as the low-temperature linear 
term in the specific heat. Agreement is found with recent measurements of 
the temperature behavior of both specific heat and pressure. These results 
suggest the presence of a very small fraction, at the parts-per-million (ppm) 
level, of two-level tunneling systems in solid 4 He, irrespective of the existence 
of supersolidity. 
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1 Introduction 

The anomalous frequency and dissipation behavior seen in torsion oscilla- 
tors (TO) [T] at low temperatures has stimulated numerous investigations, 
since it was suggested to be the signature of supersolidity [2,3,4,5,6. Re- 
cent successive TO experiments [?]I51I51Hu ] ITT P ^] confirmed the finding of the 
anomalous behavior. In addition, hysteresis behavior and long equilibration 
times have been observed [§l ll2l[T5] , which depend strongly on growth history 
and annealing i7 a . In the same temperature range, transport experiments in- 
cluding shear modulus [l"4"l[L3] . ultrasonic [TolfTT] and heat propagation [TB] 
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have shown various anomalous behaviors. However, no clear sign in mass flow 
[TMTMraiTM2MgTll221l25] and structural measurements [21125] emerged that 
can demonstrate the occurrence of a phase transition. 

It has been anticipated that thermodynamic measurements will resolve 
the existing controversy, since any true phase transition should be accompa- 
nied by a thermodynamic signature. The search for such signatures proved to 
be challenging by the experiments conducted so far, including measurements 
of the specific heat [251 27, 28, 29 . 30 f5T] . pressure dependence of the melting 
curve [321F33 . and pressure-temperature measurements of the solid [3ll35| . 
The main difficulties lie in measuring small signals at low temperatures in the 
presence of large backgrounds. With improving experiments at low tempera- 
tures, measurements down to 20 mK were conducted. While there is still no 
clear evidence of transition in the melting curve experiments, recent pressure 
measurements and specific heat measurements have both shown deviations 
from the expected pure Debye lattice behavior. Balatsky and coworkers [3"6l 
137] argued that the deviations occurring at low temperatures might be re- 
lated to a glass phase, where the role of two-level systems could be taken 
by tunneling dislocation loop segments. This description is also accompanied 
by an argument that the excess entropy associated with the deviation is too 
small for supersolidity to lead to detectable mechanical effects, if it is due to 
a supersolid fraction alone. 

In this paper, we expand the above glass model to describe the glass 
freezing transition occurring in the thermodynamic experiments of ultrapure 
4 Hc solid. We model the subsystem of tunneling dislocation segments with a 
compact distribution of two-level excitation spacings, see Fig.[T] We can then 
describe the behavior of a transition as well as the zero-temperature extrap- 
olation of the glassy behavior seen in thermodynamic measurements. Our 
results show that the low-temperature deviations in the measured specific 
heat can be explained by contributions from a glassy fraction of the solid. 
These results add further support to our previous interpretation of torsion 
oscillator experiments in terms of a backaction due to a glassy subsystem, 
which is possibly present in solid 4 He [3"51[3"71I59"] . 

2 Glass model for the specific heat 

We propose a thermodynamic model to describe the measured low-temperature 
specific heat. We postulate the existence of a distribution of two-level tun- 
neling (TLS) systems in solid hep 4 He. These TLS may be created through 
complex configurations of dislocation loops embedded in the crystal or strain 
during growth. The distributions of length and number of dislocation seg- 
ments depend on 3 He concentration, quenching and other growth processes. 
In this paper, we are comparing the effect of different growth processes on 
ultrapure 4 He containing at most (nominally) I ppb of 3 He impurities. At 
such low levels of impurities, we expect to see the intrinsic properties of solid 
4 He. 

We start with the expansion of the standard TLS model (Fig. [1} . In the 
standard glass model [40,41,36 for solids, the density of the TLS states is 
assumed to be constant, D(E) = const., to account for the linear temperature 
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Fig. 1 Density of states (DOS) of the two-level tunneling system . The black-dash 
line represents the DOS for the standard glass model 40,41,136], while the blue- 
solid line is the truncated DOS used in this work describing a two-level system with 
a cutoff energy. 

coefficient in the specific heat at low temperatures. However, it is well-known 
in the context of conventional glasses that a more careful analysis of specific 
heat data gives rise to a power law deviating slightly from linearity. The devi- 
ation may be attributed to phonon relaxation processes |42) or to a density of 
states (DOS) of the TLS that is not constant over a characteristic energy E c 
of level spacings |4*51l4"4"] . Here, we neglect the time dependence in the specific 
heat due to phonons scattering off from the TLS and for simplicity concen- 
trate on the cutoff dependence of the DOS at high energies. For example, 
such a cutoff could be due to the finite barrier height of double-well poten- 
tials giving rise to the TLS, because in real materials the tunneling barrier 
has an upper bound set by lattice and dislocation configurations [45 . This 
is also the reason why glassy behavior is usually only seen at sufficiently low 
temperatures. At high temperatures the thermal energy can easily overcome 
the barrier and the TLS effectively become a system of noninteracting single 
oscillators. At T < 0rj/5O, the specific heat of solid 4 He is well described by 

C(T) = C L (T) + C g (T), (1) 

where the phonon contribution to the molar specific heat is given by Cl{T) = 
B L T 3 , with coefficient B L = 12tt 4 R/50 3 d , R = 8314 mJ/(mol K) is the gas 
constant, and Od is the Debye temperature. The second term describes the 
glass contribution due to the TLS subsystem and is given by 

d r°° 

C g (T) = k B R— J dE EV g (E) f(E), (2) 

with ks being the Boltzmann constant and f(E) is the Fermi function. The 
DOS of the TLS may be modeled by the box distribution function 



V g (E) 



V g (E) = -V [1 - tanh((£ - E c )/W)] 



(3) 
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Here Vq is the zero-energy DOS per energy, E c is a characteristic cutoff 
energy, and W is the width of the truncated density of states. For E c — > oo, 
one obtains the standard hallmark result of glasses at low temperatures: 

C g (T) = B g T, (4) 

where B g — kuRT)®. As we will elaborate in the next section, the glass 
coefficient B g has an intrinsic finite value at low temperature even for the 
purest 4 He samples, independent of 3 He concentration. 



3 Results and discussion 

3.1 Specific Heat 

We compare our model calculations with experimental data by the Penn State 
group [301129] for four different growth processes in ultrapure He with at most 
1 ppb of 3 He impurities. We re-analyzed their experiments and obtained quite 
different assignments for the state of two of their samples: We find two solid- 
liquid coexistence samples with solid ratios 0.34 (SL34) and 0.31 (SL31) |46) . 
while we agree with the assignments for the samples grown with the blocked 
capillary method (BC). One was grown over 20 hours (BC20), the other one 
was grown over 4 hours (BC04). Notice that sample SL34 corresponds to the 
75% solid-liquid coexistence sample reported by the Penn State group and 
that the SL31 sample corresponds to their constant pressure sample (CP). 
The experimental data are described with three parameters: T>q, E c and the 
Debyc temperature Od- Wc first determine Od, or the lattice contribution, 
from the high-temperature data (see Fig. [5]) . The phonon contribution is then 
subtracted from C to obtain the difference SC = C — Cl- We fit SC/T with 
our specific heat formula in Eq.([2J for a glass. 

We show the difference in specific heat over temperature, SC/T, for 
four different growth processes in Fig. [31 The model describes well the low- 
temperature part for all four cases and the high-temperature part is within 
the scatter of the experimental data. In these plots we fixed the width of the 
cutoff to W = 1 fieV. With W <C E c , we have verified that there is no quali- 
tative difference when varying W within reasonable ranges. Notice that the 
shape of SC (T) depends strongly on the subtraction of the high-temperature 
phonon contribution. Although the Debye temperature used for the BC sam- 
ples is reasonable, it is not unique in determining the shape of SC(T), since 
a usual 1% uncertainty in the data of C(T) at high temperatures leads to 
quite large uncertainties in SC(T)/T. 

The physical and model parameters of the four samples grown under 
different conditions are summarized in Table 1. The Debye temperature Od of 
4 He increases linearly with decreasing molar volume V m in the range between 
21 and 16 cm 3 /mol [47,48 . This is also the pressure range in which the 
experimental data were taken. We find that Od for both BC grown samples 
agrees reasonably well with literature values of 28-29 K at molar volume 
V m = 20.46 cm 3 /mol and pressure P = 33 bar. The same consistency is 
found for the coexistence samples, after correcting for the solid-liquid ratio 
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T(K) 

Fig. 2 The specific heat C of solid 4 H e wi th nominally 1 ppb 3 He for four different 
growth processes taken from Lin et al. [30] . The experimental data (squares) are fit 
by BlT 3 at high temperatures T > 0.16 K (blue line). The black-dash line shows 
the total calculation with both lattice and glass contributions. 
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36.9 


BC20 


33 


20.46 


29.7 


3.0 


2.3 


6.9 


39.5 


BC04 


33 


20.46 


28.9 


6.5 


3.3 


21.5 


115.0 



Table 1 Physical and model parameters: Debye temperature On, zero energy TLS 
DOS T>o, cutoff energy E c , concentration of TLS 7itls and excess entropy AS. 



[49] . with the reported literature value of D ~ 25 K at P = 25 bar [471148] , 
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Fig. 3 5C/T for experiments (squares) and the modified glass model with a cutoff 
energy in the TLS DOS (blue line) for same four samples shown in Fig. [5J with 
SC — C — Cl- The error bars are obtained from Ref. [30] . 



A caveat is warranted regarding sample SL31. Since it has a slightly larger 
concentration of TLS than SL34, although it should be equally pure. Hence, 
one needs to consider that it was supposedly grown under CP conditions 
at P — 38 bar. However, it is not clear if the final pressure of the cell was 
actually 25 bar or higher, nor how strained the crystal was when the initial 
pressure of 38 bar was lost. All these experimental unknowns may change 
our assignment for the solid-liquid coexistence ratio for SL31. Thus its solid- 
to-liquid ratio of x — 0.31 should be considered a lower bound [39] . 

The glassy behavior is mainly characterized by the zero-energy DOS and 
cutoff energy of the TLS, which are both noticeably larger in BC04 than 
in the others. This may be explained by a rapid growth process creating a 
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strained crystal, which gives rise to both a larger TLS concentration and a 
smaller cutoff energy, i.e., a smaller maximum tunneling barrier height. On 
the other hand, the comparison of sample SL31 with BC20 leads us to believe 
that 3 He concentration does not appear to play an important role at the 1 ppb 
concentration level and below. The coexistence sample SL31 is supposed to 
be purer than BC20 ( 3 He atoms are more soluble in liquid than solid 4 He). 
Our analysis does not show significant differences between parameters for 
both cases. This is supported by our analysis that the TLS concentration of 
these samples ranges from 3.7 to 21.5 ppm, which is at least 1000 times larger 
than the nominal 3 He concentration. We argue that this small concentration 
of 3 He impurities, compared to the intrinsic concentration of TLS in even 
the best crystals grown by the Penn State group, explains the experimental 
finding that the sample grown with the BC method for over 4 hours and 
0.3 ppm of 3 He impurities exhibits nearly identical behavior as sample BC04 
with only 1 ppb of 3 He. We therefore conclude that the intrinsic property of 
solid 4 He was measured in the samples with 1 ppb 3 He impurities, given the 
much larger TLS concentration levels found. 



3.2 Entropy Analysis 

We also calculated the excess entropy, 



associated with the excess specific heat due to the low-temperature transition 
or crossover into a glass phase. The advantage of an entropy analysis over 
that of the specific heat lies in the robustness and simplicity of counting 
states in an equilibrium phase compared to detailed model calculations for 
the specific heat. We find consistently for specific heat experiments [28,29, 
130] that the obtained values are 5 to 6 orders of magnitude smaller compared 
to the theoretical prediction for a supersolid if the entire sample underwent 
Bose-Einstein condensation (BEC). In the limit of a non-interacting BEC 



one finds ASbec = 15/4(C(5/2)/<(3/2)) R{T/T c f/ 2 ~ (5/4) R ~ 10.4 J/(K 



mol). This means that if AS is indeed due to supersolidity, then the supersolid 
volume fraction is at most 11 ppm or 0.0011% in the most disordered or 
quenched sample of the four ultrapure samples studied in this work, i.e., 
sample BC04. Such a supersolid fraction in the specific heat of bulk 4 He is 
more than 100 to 1000 times smaller than is usually reported for the non- 
classical rotational inertia fraction (NCRIF) in torsion oscillator experiments. 
This enormous discrepancy between supersolid fractions in specific heat and 
torsion oscillator experiments was already noticed in Refs. [551157] . Until to 
date, this discrepancy remains a major puzzle that is hard to reconcile within 
a purely supersolid scenario. 
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Fig. 4 (a) (P - P Q )/T 2 vs. T 2 for P ~ 33 bar. The squares represent the data 
reported by Grigor'ev et al. |34II35| The black dotted line is the model calculation 
with Z) = 4.95 x 10" 3 (l/meV), E c = 0.2meV, D = 28.6K and the same W 
as that in the specific heat calculation. The blue-thick and black-thin lines are 
predictions for BC04 and BC20, respectively. The predicted curves for BC04 and 
BC20 use the same parameters as those shown in Fig. [5] (b) (P — Pq — Pl)/T 2 vs. 
T for the same data collection as in (a). This plot clearly shows that the Grigor'ev 
et al. data is characterized by a much larger Do and E c (roughly seven times larger) 
than those in the Lin et al. data. 



3.3 Comparison with Pressure Measurement 

Next we compare the specific heat with the pressure measurements and their 
deviations from a perfect Debye lattice behavior. We are using thermody- 
namic relations between the specific heat and pressure. The quantities to 
characterize the pressure measurement in the combined lattice and glass 
models are cll and a g defined by 

P(T) =P Q + P L (T) + P g {T) = P Q + a L T 4 + a g T 2 , (6) 

where P(T) is the pressure at temperature T. Po, Pl, P g are the correspond- 
ing pressure contributions of the ions at zero temperature, lattice vibrations, 
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and two-level excitations of the glass. The Mie-Griineisen theory gives the 
thermodynamic relation between pressure and specific heat 

= \t C9 > v + v LCL > v > ( 7 ) 

u / V m m 

where ji are the Gruneisen coefficients of the glass excitations (g) and of 
the lattice vibrations (L). Literature values for the Gruneisen coefficient of 
phonons in solid hep 4 He range between 2.6 < 7l < 3.0 [35. 50 , while nothing 
is known about j g of glassy 4 He. For simplicity we assume in our calculation 
of P(T) based on our model of the specific heat that 7 g ~ jl = 2.6. 

In Fig. |4]Ja) we show the temperature dependence of the pressure data 
(P — Po/T 2 ) reported by Grigor'ev et al. 04 :j 35 l as well as those predicted 
for samples BC04 and BC20 from specific heat measurements using a glass 
description. The glass contribution utilized in modeling the pressure data is 
given by 



P a {T) = 1 9 
T 2 V m 



k B R / dE EV g (E)f(E), (8) 
Jo 



with the TLS density of states D(E) defined in Eq. ©. All curves show finite 
intercepts with the ordinate, which we attribute to the glassy contribution 
on top of the T 2 -term due to lattice vibrations. The glassy contribution can 
be shown most clearly by plotting (P - P — P L )/T 2 vs. T (see Fig. @] (b)). 
The glass contribution is larger in Grigor'ev's sample and survives to higher 
temperatures compared to predictions for samples BC04 and BC20. The pa- 
rameters Do and E c are roughly seven times bigger in Grigor'ev's pressure 
data than predicted from the specific heat data by Lin et al. This may be 
expected for a much more rapidly grown solid. It is worth mentioning that 
the Gruneisen coefficient j g affects the result in the same way as Do- The 
coefficient j g can actually vary significantly from one sample to another[52 . 
Here, for simplicity, it is taken to be the same both in pressure and in specific 
heat measurements. To summarize, we have shown that the pressure mea- 
surement as well as the specific heat measurement are consistently described 
in the framework of the modified glass model. The higher TLS concentration 
extracted from Grigor'ev's data compared to Lin's data is likely due to a 
much faster cooling rate or due to a larger Griieneisen coefficient or both. 



4 Conclusions 

We have shown that a simple TLS model explains quantitatively the linear 
temperature coefficient in the specific heat and the low-temperature devia- 
tion of the specific heat, 8C, from that of a perfect Debye crystal without 
postulating the existence of supersolidity. By truncating the TLS DOS above 
a characteristic cutoff energy E c , the bump-like feature in SC can fairly well 
be captured within the scatter and uncertainty of the data. Furthermore, our 
thermodynamic analysis results in Debye temperatures 0d and TLS DOS 
T>{) that are in qualitative agreement with known P(T) measurements in the 
solid [34ll35j . This suggests the possibility of a glassy subsystem at the ppm 
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level in hep 4 He crystals. The presence of a glassy subsystem is consistent 
with recent reports of long relaxation times in torsion oscillators [H] , P vs. 
T measurements [M1I35] , and transport measurements [TOUT?] . In order to 
uniquely determine the ground state of solid He, i.e., whether it exhibits 
a supersolid or glass transition, more accurate measurements of the specific 
heat and thermal conductivity at lower temperatures and up to 0.5 K are 
needed. Finally, relaxation time measurements of heat pulses could provide 
useful and much needed information about the dynamics of a possible glass 
or supersolid subsystem in hep He solid. 
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